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ABSTRACT 

We test two different methods of using near-infrared extinction to estimate distances to dark clouds 
in the first quadrant of the Galaxy using large near infrared (2MASS and UKIDSS) surveys. VLBI 
parallax measurements of masers around massive young stars provide the most direct and bias-free 
measurement of the distance to these dark clouds. We compare the extinction distance estimates to 
these maser parallax distances. We also compare these distances to kinematic distances, including 
recent re-calibrations of the Galactic rotation curve. The extinction distance methods agree with the 
maser parallax distances (within the errors) between 66% and 100% of the time (depending on method 
and input survey) and between 85% and 100% of the time outside of the crowded Galactic center. 
Although the sample size is small, extinction distance methods reproduce maser parallax distances 
better than kinematic distances; furthermore, extinction distance methods do not suffer from the 
kinematic distance ambiguity. This validation gives us confidence that these extinction methods may 
be extended to additional dark clouds where maser parallaxes are not available. 
Subject headings: masers, ISM: dust, extinction, ISM: clouds. Galaxy: structure 



1. INTRODUCTION 

The Sun's location in the disk of the Milky Way Galaxy 
provides, in principle, an excellent vantage point to study 
the inner workings of a spiral galaxy. Stellar popula- 
tions can be resolved, chemical abundances and radial 
velocities can be obtained on a star-by-star basis, and 
the vertical, radial and azimuthal distribution of both 
stars and different phases of interstellar gas can be es- 
timated. One of the principal complications, however, 
lies in the determination of distances of individual stars, 
star-forming regions and gas clouds. The possibility of 
substantial systematic errors in distance determination 
methods have made "maps" of the Galactic stellar and 
interstellar distributions, particularly the maps of spiral 
structure, potentially unreliable. 

Complicating this process is the fact that different 
components of the Galaxy might reasonably be assumed 
to have different, albeit related, morphologies, as is ob- 
served in other spiral galaxies. |Binney fc Tremaine] 
([2008 J , for example, distinguish between bright- star 
(star-forming) arms, dust arms, mass (stellar density) 
arms, and (gravitational) potential arms, each of which 
may have a different distribution. (To these four trac- 
ers, one might add neutral gas arms, as concentrations 
of neutral hydrogen gas can have a distribution that is 
distinct from the tracers mentioned above.) 

In this paper, we concentrate on the distance to (in- 
frared) dark clouds, dense clouds of molecular gas and 
dust that are frequently associated with the early stages 
of star formation. Until recently, only kinematic dis- 
tances were available for these objects, unless they could 
be clearly associated with young stellar clusters in which 
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case photometric distances were also potentially avail- 
able. The past decade has seen the development of new 
techniques to estimate the distances to dark clouds us- 
ing the extinction of stell ar sources along the line of sight 
(e.g. M arshah et al.,2009{lKainulainen et al.j2011[ ). There 
are many ditterent possible ways of estimating extinction 
distances, with different underlying assumptions and lim- 
itations, but in general these extinction distances hold 
the promise of producing three dimensional maps of dust 
distributions independent of kinematics. 

A key point to note is that kinematic, photometric, 
and extinction distances are not absolutely calibrated. 
Although the systematic application of one method pre- 
sumably yields useful information about relative dis- 
tances, the lack of certainty about the accuracy of in- 
dividual distances has led to some justifiable skepticism 
about maps derived fro m these techniques. For example, 
Marshall et al.| ( |2009 ) found good agreement between 



kinematic and extmction distances in the first quadrant, 
but a significant (1.5 kpc) systematic difference between 
these two methods in the fourth quadrant. 

The on-going campaign to obtain radio trigonometric 
parallaxes to a set of masers across the Galactic disk 
holds the promise to test, and perhaps even to calibrate, 
other methods for estimating distance. Although the 
~ 200 parallaxes expected will not be enough to com- 
pletely map the Galaxy, they may be enough to calibrate 
these oth er methods, or at t he very least assess their re- 
liability. 



Reid et al. 



(2009), for example, have already 



developed a revised kmematic distance model based on 
early results. If the masers located in the coherent gas 
structures found in the longitude- velocity diagrams, (e.g. 
iDame et al.'2001 ) , show systematic offsets between paral- 
lax and kinematic distances, it may be possible to reliably 
map the entire gaseous structure. This would be similar 
to how photometric distances — and more recently, par- 
allax distances — of Perseus Arm star formation regions 
have been used to estimate the kinematic distance correc- 
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tion necessary to correctly map this structure (e.g. Foster 
fc MacWilliams| |2QQ6). It is important to note, tnough, 
that the intrinsic velocity dispersion of clouds will al- 
ways produce random errors in the kinematic distances 
that depend on longitude and radial velocity. Moreover, 
the near-far kinematic distance ambiguity will remain, 
and even corrected kinematic distance will be unreliable 
for sources close to / = 0° and / = 180°. 

Extinction distances have a different set of selection 
biases and systematic errors from those associated with 
kinematic distances. In this work, we use maser paral- 
laxes to test extinction distances. This is possible be- 
cause masers are generally associated with star-forming 
regions. Methanol masers, in particular, are associated 
with high-mass star formation in dark clouds. In this 
work, we present a pilot study, examining the distance 
predictions for eleven clouds between / = 5 — 60° for 
which we have near- and mid-infrared data, molecular 
line emission, and sub-millimeter observations. In ^ we 
describe our sample of clouds and their parallax distance 
measurements. The kinematic distances to these clouds 
(and their uncertainties) are discussed in ^ In Q and 
^we present two extinction distance techniques which 
we compare to the maser parallax distances. The results 
of the comparison are outlined in ^ 

2. MASER PARALLAXES 

Trigonometric parallax distances are largely free of as- 
trophysical assumptions and therefore constitute an ideal 
source of information about the structure of our galaxy. 
Recently, several groups have used Very-Long-Baseline 
Interferometry (VLBI) phase-referencing techniques to 
measure the parallax to masers in regions of massive star- 
formation with an accuracy of better than 10 //as. In the 
sample of maser parallax distance examined in this study 
the uncertainty on the distance is less than 25% out to 
6 kpc. The groups involved include the Bar and Spi- 



ral St ructure Legacy Survey (BeSSeL; [Brunthaler et al. 
2011) using the VLB A (Very Long Baseline Arr ay), a 
survey using the European VLBI Network (EVN) (jRygl 



et al. 2010), and the VLBI Explorat ion of Radio As- 
trononiyl V ERA) project (e.g. Choi et al.|2008 
eFal^[2QTl| ). 



Niinuma 



I'hese distances have been used to refine Galactic pa- 
ramet ers and rotation curves. In particular, |Reid et aL] 
([2009) used maser parallax distances to derive an up- 
dated rotation curve and found that the regions host- 
ing massive star formation are systematically rotating 
15 km s~^ slower than the rest of the galaxy. This result 
assumes that the determination of the solar peculiar mo- 
tion (motion with respect to the l ocal standard of rest 
(LSR)) is well known from [Dehnen & Binney (1998). In 
the wake of this paper, several studies have attempted 
new estimates of the solar motion, witho ut convergence 
see compa rison tables in [Francis fc An- 



to a single value 



derson 2009| Bochanski et al .^'2011'; 'Coskunoglu et al 
201 ip . However, some results (McMihan fc Binney 2010f 
suggest that the solar peculiar motion in the V direction 
(V©) was under-estimated by ~ 5 km s~^. This change 
would reduce the siz e of the lag in high -mass star-forming 
regions identified by Reid et al. (2009) to 10 km s~^, but 



not reduce it entirely. 

For this study, we select a subset of the dark clouds 
with maser parallax distances. Our methods for es- 



timating extinction distances rely upon accurately es- 
tablishing the boundaries of the densest portion of the 
cloud. Because a wealth of data exists in the 1st quad- 
rant of the Galaxy, we limit ourselves to masers in this 
region. Our ideal candidates are clouds with UKIDSS 
(UKIRT In frared Deep Sky S urvey) Galactic Plane Sur- 
vey (GPS; [Lucas et ah] 2008|), Bolocam Galactic Plane 
Survey (BGF S; [Aguirre et aL||2011) , and Galactic Ring 
Survey (GRS;p ackson et al.|2006p aata. BGPS data can 
be used to identify the boundaries of a cloud; GRS data 
can be used to assign a velocity, and to check that a 
BGPS-identified cloud is coherent in velocity; UKIDSS- 
GPS data offers higher precision near-infrared ph otome- 
try than the Tw o Micron Ah Sky Survey (2MASS; [Skrut^ 
skie et al.||2QQ6l ) . Table [l] and Table |2| list all the clouds 
and masers we consider in this stuay, although some 
methods can only be applied to a subset of these clouds. 
Additional maser parallax distan ces are continually be- 
ing reported in the literature (e. g. Nagayama et al.|2011[ 
Sanna et al. 2012 Xu et al.|2011 ) and could provide extra 



sources lor testing distance determinations. 

3. KINEMATIC DISTANCES 

Kinematic distances, estimated by assuming gas is in 
circular rotation, are simple to calculate, but have several 
sources of systematic error. Four difficulties with this 
method are (1) obtaining an accurate rotation curve, (2) 
correcting for solar motion, (3) breaking the "near-far" 
distance degeneracy in the inner galaxy (where one veloc- 
ity can map to two distances), and (4) allowing for non- 
circular motions in the vicinity of well-organized galactic 
structures. In addition, the uncertainties associated with 
kinematic distances depend strongly on longitude. 

In order to capture some of the systematic uncertain- 
ties associated with choice of rotation curve and solar 
motion, we calculate kinematic distances to the maser 
clouds using two different Galactic rota tion curves. Our 
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first set of kinematic distances use the [Clemens 
rotation curve, which includes a correction to V lsr of 
the form (7 sin/) km s~^ (which is in the same direction 
as the recent increase in of 5 km s~^). Our second set 
of kinematic distances uses the rotation curve proposed 
by iReid et all (|2009|), including the 15 km s ^ lag for 
high-mass star-forming regions (effectively a decrease of 
15 km s""*^ in Oq), since our study uses high-mass star- 
forming regions. The velocity used for these kinematic 
distances come from a variety of dense gas tracers in a 
variety of studies; see Table |2] for the velocities we have 
used. 

We estimate the random error on kinematic distances 
due to non-circular motion by perturbing the reference 
velocity by ± 3 km s~^ and deriving the kinematic dis- 
tance for the perturbed velocities. A velocity spread of 
3 km s~^ was chosen because this is the azimuthal (i.e. 
I D) velocity disp ersion of molecular clouds as reported 



Clemens (1985). Because the far kinematic distance is 



much larger than the maser parallax distance for all these 
clouds, we always adopt the near kinematic distance. We 
note that in the absence of this prior information, HI self 
absorption can sometimes be used t o brea k the distance 
ambiguity (e.g. [Roman-Duval et al.||2009 ). 



4. BLUE NUMBER COUNT EXTINCTION DISTANCES 



Maser Extinction Distance Comparison 



TABLE 1 
Masers Used in this Study 



Name 



Distance [kpc] 



Reference 



Type 



G5.89-0.39 


5.8842 


-0.3924 


1.28 


G9.62+0.20 


9.6211 


0.1959 


5.2" 


G23.01-0.41 


23.0096 


-0.4105 


4.59 


G23.44-0.18 


23.4398 


-0.1822 


5.88 


G23.66-0.13 


23.6566 


-0.1272 


3.19 


G34.39+0.22 


34.3940 


0.2215 


1.56 


G35.20-0.74 


35.1970 


-0.7431 


2.19 


G35.20-1.74 


35.2002 


-1.7364 


3.27 


W51 Main/S 


49.4884 


-0.3879 


5.41 


IRAS 19213+1723 


52.1005 


1.0429 


3.98 


G59.7+0.1 


59.7828 


0.0647 


2.16 



+ 0.09 
-0.08 
+0.6 
-0.6 
+ 0.38 
-0.33 
+ 1.37 
-0.93 
+ 0.46 
-0.35 
+ 0.12 
-0.11 
+ 0.24 
-0.20 
+ 0.56 
-0.42 
+ 0.21 
-0.20 
+ 0.67 
-0.50 
+0.10 
-0.09 



i Motogi et~aI1(|201l t 

Sanna et al. (2009) 
Brunthaler et al. (2009) 
Brunthaler et al. (2009) 
Bartkiewicz et al. ( 2008 ) 
Kurayama et al. (2011) 
Zhang et al. (2009 ) 
Zhang et al. (2009J 
Sato et al. (20^) 



Oh et al. 


(|2010 


Xu et al. 


J2OO9 



22 GHz H2O 

12 GHz methanol 

12 GHz methanol 

12 GHz methanol 

6.7 GHz methanol 

22 GHz H2O 

12 GHz methanol 

12 GHz methanol 

22 GHz H2O 

22 GHz H2O 

12 GHz methanol 



TABLE 2 

Clouds Containing Masers Used in this Study 







Cloud Velocity 








Name 


Description^ 


Ylsr [km s 1] 


Reference 


Tracer 


GRS 


BGPS 


G05.89-0.39 


UCHII 


+7 


Arikawa et al."('1999") 


CO 


N 


Y 


G09.62+0.20 


HCHII 


+4 


ScovUle et al. ( 1987 ) 


CO 


N 


Y 


G23.01-0.41 


HMSFR 


+77 


Codella et al. (1997) 


NH3 


Y 


Y 


G23.44-0.18 


HMSFR 


+99 


Codella et al. (I^UIU) 


NH3 


Y 


Y 


G23.66-0.13 


HII region? 


+80 


Urquhart et al. (l^Ull ) 


NH3 


Y 


Y 


G34.39+0.22 


HMSFR 


+57 


"Bfontman et al. ( 1996 ) 


cs 


Y 


Y 


G35.20-0.74 


HMSFR 


+35 


Solomon et al. ( 19«7 ) 


CO 


Y 


N 


G35.20-1.74 


UCHII 


+43 


STutzki et al. (19«4) 


NH3 


N 


N 


W51 Main/S 


HMSFR 


+55 


Zhang et al. ( 199« ) 


cs 


Y 


Y 


IRAS 19213+1723 


HII region? 


+42 


Hrontman et al. ( 1996 ) 


cs 


Y 


N 


G59.78+0.06 


HMSFR 


+29 


"Bronfman et al. ( 199t) ) 


cs 


N 


Y 



^ UCHII Ultra-Compact HII Region, HCHII Hyper-Compact HII Region, HMSFR High-Mass Star-Forming 
Region. 



Our first extinction distance method, the Blue Number 
Count method, uses the number of foreground stars and 
a Galactic model of the stellar distribution to estimate 
the distance to dark clouds. Using star counts to esti- 
m ate t he distance to dark clouds goe s back at least as far 
aslmifK 1923) and "Bok" (1931) (see Krelo wski fc Papaj 
Il993[ tor a history of the origins of this field). Our version 
oT this method exploits the relatively narrow range of in- 
trinsic near-infrared colors of stars to separate foreground 
from background stars. Clouds with sufficiently high col- 
umn density will redden background stars to such an ex- 
tent that the intrinsically reddest foreground star will be 
bluer than the intrinsically bluest background star. We 
can then simply count the number of stars bluer than 
this color threshold and compare with a Galactic model 
of the stellar distribution to estimate a distance. This 
simple method has the advantage that the distribution 
of stellar colors in the Galactic model need not be pre- 
cisely correct, only the number density as a function of 
distance from the Sun. This method is quite similar to 
the method applied to derive dist ances to n e arby (D < 
1 kpc) molecular clouds in Lom bardi et al.| (2Q1Q|) and 
to more distant clouds in ^Kainulainen et al.| (|2Ull|) for 



distant Infrared Dark Clouds (iKDCs). 

The Blue Number Count method involves three steps: 
(1) define the region of the dark cloud associated with a 
maser source which is dense enough produce a clear color 



separation; (2) count the number of blue stars within this 
region, accounting for confusion and incompleteness; (3) 
compare the number of blue stars with a Galactic model 
to determine a distance to the cloud. In the following 
sections we elaborate on each of these steps and illustrate 
the method in Figure [T] |2j and |3] for an example cloud, 
G23.01-0.41. We also discuss some of the possible sources 
of bias and error in this method. 



4.1. Defining a Dense Region 

The first step in the Blue Number Count method is to 
identify the column density necessary to make a clear 
separation between foreground and background stars. 
The range of int rinsic J — K colors is less than 1.5 
magnitudes (e.g. [Lombardi fc Alves 2001|. Therefore 
we need enough column density to produce at least 1.5 
magnitud e s of r eddening m. J — K. Using the Rieke fc 
Lebofsky (1985) extinction curve Aj/Ay = 0.282 and 
Ak/Av =TnT2, thus Av = A{J-K)/0.170 and the ex- 
tinction required to produce the requisite amount of red- 
dening to clearly separate foreground from background 
stars is 8.8 magnitudes of Ay. 
We use Bolocam Galactic Plane Survey (BGPS; 



Aguirre et al. 2011) data to identify the boundaries of 
the cloud with column densities large enough to produce 
an extinction of at least 8.8 magnitudes of Ay. We as- 
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Fig. 1. — UKIDSS three-color image (blue = J, green = if, 
red=K) of G23.01— 0.41. The yellow contour shows the bound- 
ary of this cloud as determined from the BGPS thermal emission 
map. The cyan cross marks the position of the maser which has 
been used to measure a parallax distance toward this cloud. Note 
that inside the contour the stellar colors are typically blue or very 
red^s confirmed by the J — K color histogram for this cloud (Fig- 
ure [2}. 
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Fig. 2. — J — K color histogram for stars within the contour 
defined for G23. 01-0.41 using UKIDSS data. This cloud produces 
a very clear separation between blue foreground stars (J — K < 
1.5 mags; dark gray) and red background stars (J — K > 1.5 mags; 
light gray). The results are insensitive to the exact value of the 
color cut as foreground stars are well separated from background 
stars. 

sume that 



I-beam 



H2 



(1) 



where S^^^^ is the flux per beam, is the beam sohd 
angle, /iHs is the mass per hydrogen molecule, hz^ is the 
dust opacity, and By{T) is the Planck function. For 
Bolocam data, A = 1120 /im, and Ohpbw = 31''. We 
used /iH2 = 2.8, and f<ijy = 0.0114 cm^ correspond- 
ing to the wavelength-interpolated fOssenkopf fc Henning 
( |l994J value for dust grains with thin ice mantels in dense 



1000 2000 3000 4000 5000 6000 7000 8000 
Distance [pc] 

Fig. 3. — Blue Number Count distance to G23.01-0.41 using 
UKIDSS data. We compare the (completeness-corrected) number 
of blue (J — K < 1.5 mags) stars per acrminute within the cloud 
contour (horizontal line) against the number of stars predicted to 
be closer than a given distance from the Besangon Galactic model 
of Robin et al. (2003) (dashed line). Poisson errors on the number 
of stars produce an estimate of the (non-systematic) uncertainty 
of this distance determination (vertical lines). 

regions (n = 10^ cm~^) and assumed a dust temperature 
of T = 10 K. For these assumptions 

NhJSI''^'^ = 6.77 X lO^^cm-^py beam-^)-^ (2) 

in the appropriate units (since BGPS data is available in 
maps with units of Jy beam~^). N(H2) is then converted 
to Ay assuming 



= 9.4 X 10^^ cm- \Av/ mag) 



(3) 

from "Bohlin et al.l (^978^. Therefore, our critical flux 
threshold ol 8.8 magnitudes corresponds to a BGPS flux 
threshold of 0.12 mJy beam~^. The various conversion 
factors assumed here (T, a^j^, and the relation between 
Ay and N(i^2)) are relatively uncertain, but the exact 
value we obtain is not critical to our analysis; all we re- 
quires is that the adopted flux threshold produce a sep- 
aration between foreground and background stars. Fig- 
ure [l] displays this 0.12 mJy beam~^ contour for our ex- 
ample cloud (G23.01— 0.41) and shows how this contour 
matches the area over which stars are either very blue (in 
front of the cloud) or very red (behind a large column of 
dust). 

Three of our regions are not included in BGPS because 
their Galactic latitudes were too large to be included in 
the BGPS. For two of these sources (IRAS 19213+1723, 
G35.20-0.74) we use the GRS ^^CO data, integrated over 
±15 km s~^ and convert this to a column density using 
the conversion factors in Carpenter (2000). For G35.20- 
1.74, no GRS or BGPS data is available, and we estimate 
the Ay =8. 8 mag region by hand, tracing the obvious 
dark extinction area. A few other cloud boundaries are 
adjusted to exclude regions of particularly bright emis- 
sion or obvious foreground clusters of blue stars (see Fig- 
ure |4|. 

4.2. Counting Blue Stars 



Maser Extinction Distance Comparison 
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Fig. 4. Contour for G35.20-0.74 made from the GRS ^^CO 
integrated intensity image overlain on the 2MASS three color image 
(J=blue, if=green, K=red). The contour was adjusted by hand 
to exclude the bright blue cluster of stars at ^=35.12°, fe=-0.75°. 
This cluster may have been born in this cloud, and its exclusion is 
necessary so that the distance determination is not biased. Note 
that embedded young stars such as those around the maser (cyan 
cross) will not bias the Blue Number Count method, since red stars 
are not used in the analysis. 

The second step in the Blue Number Count method 
is to count the number of blue stars within the region 
of interest while correctly accounting for completeness. 
The masers reside in regions of high extinction and ac- 
tive star formation. Active star formation can produce 
significant nebulosity and thus a spatially variable back- 
ground. On the other hand, the high extinction means 
that confusion is less severe when performing photome- 
try on a dark cloud. Average completeness estimates are 
therefore unreliable; the completeness of a given survey 
is different for each cloud. 

We performed our own estimates of the completeness 
by injecting synthetic stars at random positions within 
the Ay =8.8 mag cloud contour and testing for their re- 
cove ry usin^ automatic phot ometry with Source Extrac- 



tor (Bertin & Arnouts |1996 ). For example, to estimate 
the completeness ot a UKiDSS field we set the satura- 
tion level at 40,000 ADU and used a 5-pixel (2^^) radius 
photometric aperture. This size aperture corresponds to 
the default aperture size APERMAG3 in the UKIDSS cat- 
alog. We ran 500 trials on each cloud, injecting stars 
as artificial circular gaussians into the i^-band UKIDSS 
image. We estimate the median FWHM from all the 
stars in the image and use this value for the FWHM for 
our artificial stars. The noise comes from sitting on top 
of the noisy background of the image. In each trial we 
inserted stars with uniform integer magnitudes between 
11 and 19 magnitudes, inclusive. We tested whether we 
retrieved input stars by searching for output stars within 
1 pixel of the input position with a magnitude difference 
of no more than 0.5 magnitudes plus the magnitude un- 
certainty reported by Source Extractor. In addition we 
required a Source Extractor FLAG value of less than 4, 
eliminating stars with saturated pixels or truncated ob- 
jects but allowing the star to be partially blended with 
another star. 

We estimated the completeness of the field at integer 



magnitudes and interpolated between these values in or- 
der to weight each measured star appropriately, as shown 
in Figure [5| Because the completeness approaches zero 
beyond 17 magnitudes in we limit the stars consid- 
ered to be between of 11 and 17 magnitudes in order to 
avoid applying very large correction factor. We estimate 
completeness from the i^-band image since the method 
only uses the blue stars {J — K < 1.5 mags), and the 
J image is complete to significantly fainter magnitudes. 
A similar procedure was carried out for 2MASS images; 
for 2MASS we limited our analysis to 7 mag < K < 13 
mag. Our estimate of the completeness in the 2MASS 
images is brighter than the quoted survey completeness 
for uncrowded refflon s of 2MASS (14.3 magnitudes in K; 
Skrutskie et al.||2006[ ) because (1) even with dark clouds 
blocking out niany background stars there are enough 
foreground stars in the Galactic plane to cause some con- 
fusion, (2) there is significant near-infrared nebulosity 
in many regions, and (3) our parameters for identifying 
sources and performing photometry may be less sensitive 
and tolerant of source-blending than the procedures used 
to generate the 2MASS catalog. 

We can not use the standard point-source catalogs from 
UKIDSS or 2MASS because our completeness estimate 
is only appropriate for our particular photometric pa- 
rameters. Instead, we create our own catalog with the 
same Source Extractor parameters and quality checks as 
used in the completeness study and fit these measured 
magnitudes against UKIDSS and 2MASS magnitudes to 
derive magnitude zero-points which puts our photom- 
etry roughly on the UKIDSS and 2MASS photometric 
systems. Because our results do not rely on exact col- 
ors, small color or magnitude offsets with respect to the 
standard catalog magnitudes are unimportant. Figure [2] 
shows the J — K stellar color histogram for our exam- 
ple cloud (G23.01— 0.41). The color separation between 
foreground and background stars is quite pronounced. 

4.3. Comparing to a Galactic Model 

The final step in the Blue Number Count method is to 
compare our count of foreground stars with a Galactic 
stellar distribution model. We use the Besancon model 
dRobin et aL]|2003 ) to generate a synthetic stellar cata- 
log consist mg of distances and near-infrared magnitudes 
for each of our target fields. The Besancon model al- 
lows the generation of a "small-field" catalog of arbitrary 
size. That is, we can obtain an arbitrarily large sample 
of stars at the desired position to reduce Poisson noise in 
the simulated catalog. We assume that gradients in the 
stellar distribution across our fields are negligible. We in- 
clude all spectral types and populations available in the 
model and initially include no diffuse extinction. The 
model output is a catalog of stars containing both in- 
trinsic properties (i.e. distance) and observed quantities 
(i.e. near-infrared magnitudes). 

We process the catalog produced from the Besancon 
model to include diffuse Galactic extinction. We use the 
default model value of 0.7 magnitudes of Ay per kpc. 
Since typical maser distances are between 1 and 5 kpc, 
this gives us less than 3.5 magnitudes of visual extinction 
from the diffuse component which translates into a max- 
imum E{J — K) of 0.57 magnitudes. Our method does 
not work if there is a foreground cloud which produces 
significant excess extinction and our method should only 
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Fig. 5. — X-band completeness estimates for our sample of clouds containing masers for [left] UKIDSS and [right] 2MASS. Completeness 
was estimated by injection and recovery of 500 artificial stars at each integer magnitude. All artificial stars were placed randomly within 
the cloud contour. Due to the steep drop-off in completeness, we limited our analysis to stars brighter than 17 magnitudes in UKIDSS and 
13 magnitudes in 2MASS. 



be applied to clouds where auxiliary information con- 
firms this exclusion. Small changes to the diffuse ex- 
tinction with distance will produce a small change in 
the distance estimate for nearby clouds and a moderate 
change for distant clouds. If the true diffuse extinction is 
larger than we assume, then some foreground stars could 
be reddened so much that we incorrectly identify them 
as background stars. The exclusion of these foreground 
stars would make clouds appear closer than they really 
are; this is of particular concern for distant clouds where 
the impact of diffuse extinction is more significant. 
An alternative approach to dealing with diffuse extinc- 



tion is adopted by Kainulai nen^t a l.| (|2Qn )^ who use the 
3D extinction profiles of Marshall et ctE " pt)Q6| ) to model 
the extinctio n as a function of dist ance along each line 
of sight. The [Marshall et 




using the same 
near-infrared data (2MASS) 
to disentangle the uncertainty in these profiles from the 
other uncertainties in the model. By focusing only on the 
very dense portions of clouds we are able to reduce our 
sensitivity to the exact value of the diffuse extinction at 
the expense of using a smaller region of the cloud for our 
analysis. As long as there is a distinct break in the col- 
ors of stars we can separate foreground from background 
stars, and our only dependence on the Galactic model is 
the number density of stars as a function of distance. 

We use the model catalog to define a monotonically 
increasing relationship between the cloud distance and 
surface number density of foreground (blue: J — K < 
1.5 mags) stars per square arcminute. This relationship 
saturates at some distance due to a combination of dif- 
fuse extinction and (depending on Galactic position) the 
finite size of the Galactic disk, but remains monotonic 
for all our clouds well past the estimated distance. We 
solve for the distance to a cloud by finding the distance 
at which our synthetic curve of blue foreground stars 
matches our estimated true number of blue stars, after 
correcting for completeness. We assume that the major 
source of random (non-systematic) uncertainty is Poisson 
{^/n) noise due to our small number of blue foreground 
stars and use the number of measured blue stars (before 
correcting for completeness) in this calculation. Figure [s] 



shows this distance determination for our sample cloud 
(G23.01-0.41) using UKIDSS. 

5. RED GIANT EXTINCTION DISTANCES 

Our second method is the Red Giant extinction dis- 
tance method. This method makes extinction measure- 
ments to individual giant stars along the same line of 
sight as a molecular cloud. If there are enough fore- 
ground and background giant stars available then the 
presence of the cloud will be apparent in an extinc- 
tion measurement histogram. An appropriate extinction- 
distance relationship is then used to convert these mea- 
sured extinction values into a distance. One of the 
strengths of this method is that if enough field stars 
are available then it is possible to detect multiple clouds 
along the same line of sight. 

A fu ll descrip tion o f this method was presented in 
Stead fc Hoare] ( 2Q1Q[ ). However, the most important 
aspects ot the method are also detailed below for com- 
pleteness. The photometric error cuts for each survey 
are discussed below. However the UKIDSS data are also 
subject to further quality control cuts to remove stars 
with pixels close to saturation, blended objects and other 
artifacts that may affect the photometry. For further de- 
tails on the UKIDS S data quality controls please refer to 
Lucas et all (120081). 



5.1. Field Star Selection 

The size of the region of sky used to extract field stars 
depends upon several factors. The area of sky must be 
large enough to yield an adequate number of extinction 
measurements to both foreground and background stars. 
However, the larger the area used, the greater the con- 
tamination in the extinction histogram. This contami- 
nation occurs predominantly due to the column density 
structure of the molecular cloud, however overlapping 
clouds and other irregularities in the interstellar medium 
will also contribute. 

To minimize the contamination in the extinction mea- 
surements, where possible we have used GRS data to de- 
fine the area of sky used for extraction. Rat hborne et al. 
(2009) identify molecular clouds in the GRS, specifying 
both their spatial position and their Vlsr- Integrating 
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Fig. 6. — The Red Giant extinction method for G23.01— 0.41. [Left; a]: an Ay histogram generated using UKIDSS data extracted from a 
region containing the maser G23.01-0.41. The near and far sides of the molecular cloud are marked with dashed lines, derived from the la 
deviations of the skewed Gaussians fitted to each histogram (solid curve). [Center; b]: same as (a) except using 2MASS data. [Rig ht; c1_ : 



A distance-Ay plot showing the extinction-distance relationship (solid thin line with error bars) created using the Marsha ll et al.| (I2006| 
data centered at G23.0-0.50. The two sets of dashed lines illustrate how the near and far sides of the molecular cloud, the dasJied Imes 
from (a), are converted into a distance and associated error. 



the ^^CO luminosity measured in the GRS along a spe- 
cific range of Vlsr will map out the molecular cloud 
column density. The GRS ^^CO maps give the ability 
to select an area of sky containing a similar cloud col- 
umn density. Isolating these regions will therefore limit 
the uncertainty caused by cloud patchiness. As previ- 
ously stated, the exact size of the region should be as 
large as possible while minimizing the noise in the ex- 
tinction histogram. The size of each region is decided on 
a cloud-by-cloud basis after assessment of each extinction 
histogram. 

5.2. Isolating Red Giants 

Late type giant stars are old enough to have undergone 
several Galactic rotations. They have therefore moved 
from their original birth places in the spiral arms, with 
reference to the molecular clouds, and they can be consid- 
ered ubiquitous throughout the entire Galaxy. Late type 
giants are very numerous and are also infrared bright. 
These qualities, plus the fact that they have known col- 
ors, are why many authors have used them to create ex- 
tinction ma ps of the Galaxy (Marshall et al.^^2QQ6^ Lom^ 
biIdi||2QQ9l ). 

Late type giants occupy well known regions on 
both color-color and color-mag nitude d iagrams 
(^CCDs/CMDs). 



. ^ . ^ For this reason |Stead fc Hoare| 
f2010) were able to use a series of color cuts to isolate 
them from the general field population. They tested 
their method, with realistic synthetic data, to confirm 
that over 90% of the final sample of stars will be GOIII 
stars or later. Their method has been replicated for the 
work in this paper. 

5.3. Extinction Measurements and Ay Histograms 



In the same manner as Stead & Hoare (2010), the 
color-selected sample of late type giants have been de- 
reddened, in color-color space, using a KOIII reddening 
track to their point of intersection with the giant locus. 
This reddening track has been created u sin^ a KOIII stel- 
lar model from Castelli fc Kurucz ( 20041) and th e average 

^- - L- I _i I LU I I _ -L i /foTTTTTTl _ o 1 /i\ 



extinction power law ot Istead fc 



Hoare (2009, a=2.U) 



using the method described therein. 

The length of the reddening track traversed during the 
de-reddening process relates directly to the amount of 



extinction each individual star suffers. Photometric er- 
rors, as well as directly adding to the error in the extinc- 
tion measurement, will also make the point of intersec- 
tion with the giant locus more uncertain. Photometric 
errors are the greatest source of noise in an extinction 
histogram. The faintest late-type giant stars in each re- 
gion, i.e. those with the largest photometric errors, will 
be those stars that are behind each molecular cloud. As 
such stars are essential to the detection of each molecular 
cloud there is therefore a tradeoff between completeness 
and noise in the extinction histogram. For the purpose of 
this paper we consider only stars with photometric errors 
of less than 0.05 mags in the UKIDSS catalog and less 
than 0.10 mags in the 2MASS catalog. The less strin- 
gent error limit on the 2MASS catalog is used to provide 
a sufficient number of stars. 

A value of A j is determined directly for each individual 
star using this de-reddening process. This value is then 
scaled to a value of Ay using the ratios A j/ Ay =0.2833 
and Aj/Ay=0.2899 for UKIDSS and 2MASS respec- 
tively. These have been calculated using the Cardelh] 
indl|l989) extinction curve with Ry=3.1. As noted in 



the literature, the value Ry=3.1 may not be appropriate 
for molecular clouds, and our extinction is likely pro- 
duced by a combination of diffuse extinction and extinc- 
tion from molecular clouds. Adopting an extinction law 
is necessary, since we must compare our derived values 
of Aj with the extinction distance relations of Marshah] 



et al.| ([2006) which are in A^^. However, the uncertainty 



mtroduced by our choice of Ry is small. The conversion 
between Aj and A^^ differs by only 1.3% if we choose 
Ry = 5.5 instead of Ry = 3.1. Note that since our choice 
of Ry simply determines the multiplicative conversion 
factors between the different wavelengths, the fact that 
we are converting Aj and A^^ to Ay before making 
the comparison does not introduce any additional uncer- 
tainty. 

We present an example of this method using the same 
cloud (G23.01— 0.41) as was used to illustrated the Blue 
Number Count method (Figure [6|. The cloud that con- 
tains the maser was mapped in "CO in the GRS. Stars 
have been selected from within the >8 K km s~^ GRS 
^^CO integrated intensity contour. As mentioned pre- 
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viously, this value was selected after inspection of the 
Ay histogram to optimize the completenes s /noise tr ade- 



off. Using the conversion factors in , Carp enter (20001) this 
corresponds to an Ay of 4.1 magnitudes, ana is thus a 
lower column density threshold than used in the Blue 
Number Count method. 

Figure [g] (a) displays an Ay histogram created using 
UKIDSS aata. The histogram begins with a small rise 
in star counts, peaking at Ay~3, and then quickly dips 
again to a minimum at Ay~5. Following this dip the 
star counts quickly rise again to form a second peak at 
Ay ^10. This double peaked histogram is consistent with 
the idea that there is a single molecular cloud along the 
line of sight. Stead & Hoare ( |2010 ) define the gap be- 
tween the two peaks, measured by fitting skewed Gaus- 
sians to each peak to determine the respective near and 
far 1(7 widths, as the Ay of the cloud. They used syn- 
thetic data to determine the error in this measurement 
to be Ay±0.5 mag. This value has been included in all 
of the distance error assessment that follows. 

Figure|6](b) shows the same Ay histogram for 2MASS 
data. Trie two obvious differences between each his- 
togram are the higher star counts in the deeper UKIDSS 
data and the more jagged 2MASS histogram, owing to 
the smaller number of stars and the less precise 2MASS 
photometry. 

5.4. Distance-Ay plots 

'Marshall et al." (2006) used the Besancon Galactic 
model and 2MASS data to generate a three dimensional 
map of the Galactic interstellar extinction distribution. 
They gridded the Galaxy into 15^x15' tiles and gener- 
ated an extinction-distance relationship (EDR) for each 
tile. The line of sight extinction measurements in this pa- 
per are converted to distances using the Mars hall et al . 
|2006) EDRs. The spatially closest EDR to each cloud 
has been used. 



The Marshall et al, 



( 2006 ) data are presented in terms 
of distance vs. A^^, the extinction in the 2MASS Ks 
filter. Figure J6 (c) contains the EDR centered at G23.0- 
0.50. The UKIDSS measurements produce an estimate 
for the distance to the cloud of D = 4.3^0 3 ^P^- This 
process is also repeated using 2MASS data yielding a 
distance, D = 4.3lo.3 ^P^- Typically the errors produced 
using 2MASS data are larger than those derived using 
UKIDSS data. For this cloud however the 2MASS data 
have constrained the result somewhat more reliably than 
with UKIDSS data. The reasons for this are discussed 
below. 

5.5. UKIDSS versus 2MASS 

As previously discussed, an adequate number of both 
foreground and background stars are required to make 
the Red Giant method work. For this reason if a cloud 
is too close, and there are too few foreground stars avail- 
able, it will not be possible to constrain a distance to the 
cloud. Likewise, if the photometry is not deep enough, 
too few background stars will be detected. However 
whereas the ability to derive distances using 2MASS data 
depends only upon photometric depth, the UKIDSS data 
have a completeness issue at the brighter end caused by 
the saturation limit. Stars in UKIDSS typically satu- 
rate at mag when selecting only the most reli- 
able data. Therefore, for nearby clouds most foreground 



giant stars may be saturated and thus unavailable in 
UKIDSS. UKIDSS and 2MASS complement each other 
well in this respect, with 2MASS able to constrain dis- 
tances to nearby clouds and UKIDSS able to probe to 
much further distances. In the above example, towards 
this particular region of the Galaxy, ~4 kpc is the point 
where 2MASS outperforms UKIDSS. This is apparent in 
Fig. [6] (a) and (b) when a comparison is made between 
the number of foreground stars in each histogram. A 
merged catalog containing 2MASS and UKIDSS data for 
the bright and faint stars, respectively, would be ideal, 
but this is beyond the scope of this work. 

6. COMPARISON OF METHODS IN 1ST QUADRANT 
SOURCES 

Table [3] lists the distance determinations for all meth- 
ods. We do not report a distance measurement for some 
combinations of sources and methods, either because the 
method fails to identify a cloud or because certain data 
are unavailable. For instance, we do not report a Blue 
Number Count method distance if there are fewer than 
five blue foreground stars. For clouds with small angular 
extent there are often not enough blue 2MASS stars to 
make a measurement. For G35.20— 1.74 UKIDSS data 
images were not available, although catalog information 
was available. Therefore, the Red Giant method (which 
uses the catalog) was able to determine a distance for 
this source using UKIDSS data, but the Blue Number 
Count method (which requires the image data to esti- 
mate completeness) was not able to mak e a determina- 
tion using UKIDSS. As discussed in §5.5[ the Red Giant 
method requires foreground and background red giants, 
so for some clouds too many of the foreground giants are 
saturated in UKIDSS; conversely, in some clouds there 
are insufficient background red giants detected in the less 
sensitive 2MASS catalog. 

Figure [jl and Figure [S] compare all our methods for our 
sample oflst quadrant masers sorted by Galactic longi- 
tude and maser parallax distance, respectively. The two 
sources near the Galactic center are shaded red in both 
these figures. For these sources kinematic distances are 
unreliable and there is strong evide nce for at least two 
clouds along the line of sight (see ^6.1.1 for more de- 



tails). As expected, the different distance methods often 
disagree for these two clouds. Table [4] shows the absolute 
number and fraction of sources for which each method is 
within 2(7 agreement of the maser parallax distance (us- 
ing the a for each method added in quadrature) . As this 
table shows, the methods generally agree within their 
errors. 

There are a few notable exceptions where the methods 
fail. These clouds are G23.66-0.13, G34.39+0.22, and 
the two clouds near the Galactic center (G05.89— 0.39 
and G09 .62+0.20). We discuss these clouds individually 
For the rest of the sources the agreement is 



6.1 



generally good. In particular there are clouds (such as 
G35. 20-0.74 and G23.01-0.41) where ah the distances 
agree extremely well and have small uncertainty {a<0.5 
kpc). 

Figure [8] suggests that there is a systematic tendency to 
overestimate the extinction distances for nearby sources. 
In particular, G05. 89-0.39, G34.39+0.22, and (to a 
lesser extent) G59. 78+0.06 all have maser parallax dis- 
tances < 2.5 kpc and kinematic/extinction distances 
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TABLE 3 
Distances to All Clouds 



Name 


Maser Parallax 


Blue Number Count Extinction 


Red Giant Extinction 


Kinematic 






2MASS UKIDSS 


2MASS 


UKIDSS 


Clemens (1985) 


Reid et al. (2009) 


G05.89-0.39 


-1 90+0.09 

J..^O_Q Qg 


4 5+0-6 




q 7 + 0.3 

6 2+°-^ 


Z.U_Q 7 




G09.62+0.20 






44+0.4 
7 4+0-7 








G23. 01-0.41 




4.0l?;9 4.7t°i 


4 3+°-'^ 
^"^-0.3 


4.3l°J 




4 6+0-1 


G23.44-0.18 


00 + 1.37 
O.»»_o.93 


c: 1+0.6 A 0+0.2 
^•-•--0.7 4-^-0.2 


5 6+0-3 


4.8ii:? 


^ 0+0.2 


54+0.1 


G23.66-0.13 




4 9+1-3 








4 7+0-1 
^•'-0.1 


G34.39+0.22 


-L-OO_0.il 


9 Q+0.9 q ^+0.3 


2 6+°-^ 

q q+0.6 


2 6+°-^ 

^+0.8 


q 7+0.2 
' -0.2 


q p^+0.2 


G35.20-0.74 


^•-^^-0.20 


9 .+0.8 9 7+0.3 
^•^-0.7 ^-'-0.3 


94+0.3 


2 6+°-'^ 


r> q + 0.2 

^•'^-0.2 


2 4+°-2 


G35.20-1.74 


q r)7+0.56 

•J-^' -0.42 


2 5+°-'^ 

^•^-1.0 






9 0+0.2 


9 0+0.2 


W51 Main/S 


541+0.21 
^•^-■--0.20 






5.9ii:i 




5 5+0.8 


IRAS 19213+1723 


q Qo + 0.67 






4.61^:^ 


q p:+0.4 


q 0+0.5 


G59.78+0.06 


2 16+°-^° 

^•-^"-0.09 


3.81^2 


9 p:+0.2 




q q+1-0 
•^"^-O.S 


4 2+^-° 



Note. — All distances are in kpc 



TABLE 4 

2(7 Agreement with Maser Distances 



Method 


Full Sample 


Excluding ^ < 15 


Blue Number Count (2MASS) 


6/6 (100%) 


6/6 


(100%) 


Blue Number Count (UKIDSS) 


6/9 (66.7%) 


6/7 


(85.7%) 


Red Giant (2MASS) 


7/8 (87.5%) 


6/6 


(100%) 


Red Giant (UKIDSS) 


7/9 (77.8%) 


6/7 


(85.7%) 


Kinematic (Clemens) 


7/11 (63.6%) 


6/9 


(66.7%) 


Kinematic (Reid) 


8/11 (72.3%) 


7/9 


(77.8%) 



which are systematically larger. Nearby clouds can 
present particular challenges for both extinction meth- 
ods. If all foreground giants are saturated, then the Red 
Giant method will fail to see a clean separation. Very 
nearby clouds also have a low density of foreground blue 
stars for the Blue Number Count method, although this 
is often offset by the fact that nearby clouds are larger 
in angular extent. 

Both extinction distance methods rely on average stel- 
lar distributions, so small-scale inhomogeneities in the 
stellar population can produce erroneous distance esti- 
mates. In the Blue Number Count method the reliance 
on average stellar properties is explicit since we use the 
average stellar density along a given line of sight; in the 
Red Gian t method this d ependence arises through the 
use of the Marshall et al.| ([2006) extinction-distance re- 
lationships, which rely on average stellar colors. Small- 
scale inhomogeneities will produce an incorrect distance 
in the Blue Number Count method; for instance, a (weak, 
hard-to-detect) stellar cluster projected onto the line of 
a dark cloud would cause us to estimate an incorrectly 
large distance. For more distant clouds, we expect that 
small-scale inhomogeneities in the foreground population 
will tend to average out and be less significant. The 



extinction-distance relations of Marshall et al. (2006) 
used in the Red Giant method rely on fitting the full 
color distribution of stars toward a given position. Since 
this method considers all stars toward a given line of 
sight, we would not expect small-scale inhomogeneities 
to produce a bias for any particular distance. 

One possible explanation why both the extinction and 
the kinematic distances fail for certain clouds would be 
the presence of a spiral arm. Spiral arms may be as- 
sociated with shocks and therefore departures from the 
circular Galactic rotation curve. In addition, spiral arms 
are not included in the Besancon Galactic model, and 
so neither the stellar populations in these arms, nor the 
extinction within them, are properly modeled. The old 
stellar population (including the red giants) will typically 
be well mixed, and it is only the young stellar population 
that contributes to the arm/inter arm contrast. Spiral 
arm tangencies in the first quadrant are around / = 25° 
and / = 35° ( Benjamin 2009) . Unsurprisingly, we have 
three clouds along each of these lines of sight, as regions 
of high-mass star formation tend to occur in spiral arms. 
Most of these clouds have extinction distances that agree 
with the maser parallax distances, suggesting that look- 
ing along spiral arms does not tend to produce erroneous 
distance estimates. 

Finally, it is possible that the maser parallax distances 
may be incorrect. Although the maser parallax dis- 
tances require significantly fewer astrophysical assump- 
tions, there are possible situations in which these paral- 
lax determinations could fail. The fundamental assump- 
tion in these observations is that the masers in these ob- 
jects persist and move in a regular fashion over the sev- 
eral years when observations are made. This is typically 
mitigated by the use of multiple maser spots. However, 
uncompensated atmospheric delays can produce system- 
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Fig. 7. — Comparison of different distance estimates to the maser clouds, sorted by Galactic longitude. Black lines and gray shaded 
regions show the maser parallax distance la uncertainties taken from the literature. Blue points show Blue Number Count distances using 
2MASS (open symbols) and UKIDSS (filled symbols). Red points show Red Giant distances using 2MASS (open symbols) and UKIDSS 
(filled symbols); this method sometimes identifies multiple clouds along the line of sight, in which case both distances are shown. Green 
points show kinematic distances based on two different rotation curves. All error bars are la and are estimated for each method as explained 
in the text. Red shaded clouds show the two sources near the Galactic center where kinematic distances are unreliable and there is evidence 
for multiple clouds along the line of sight. Figure [s] shows this figure sorted by maser parallax distance. 



atic errors in the parallax measurement. 

There are a few masers for which multiple groups have 
measured the parallax independently, and in general the 
agreement is within the estimated uncertainties. How- 
ever in one case, G48.61+Q. Q2, the two distance mea sure- 
ments are lOa discrepant (Nagayama et al. 20111 Reid 



private communication; ). I'his demonstrates that it is 
possible that some maser parallax distances are incorrect. 

6.1. Comments on Individual Sources 

6.1.1. Galactic Center Sources (C05.89-0.39 and 
C 09. 62+ 0.20) 

It is difficult to determine distances for sources that 
are close to the Galactic center. Kinematic distances 
are unreliable close to the Galactic center since most or- 
bital motion is perpendicular to the line of sight. Both 
extinction methods have problems if there are multiple 
clouds along the line of sight; the projected surface den- 
sity of molecular clouds, and thus the chance of multiple 
clouds along the line of sight, increases toward the Galac- 
tic center. With sufficient data, the Red Giant method 
can distinguish multiple jumps in extinction and there- 
fore identify multiple clouds. The Blue Number Count 
distance method is sensitive only to the ffist dense cloud 
along the line of sight. 

For the two sources close to the Galactic center 
(G05.89-0.39 and G09.62+0.20) one or more of the 
methods presented in this paper fails to obtain a distance 
consistent with the maser parallax distance. The Red 
Giant method identifies multiple clouds along the line 



of sight toward both these positions. Spectra from the 
Millime ter Astronomy L egacy Team 90 GHz (MALT90) 
survey ( Foster et al.|201 1) conffim that there are at least 
two clouds by detecting two velocity components in dense 
molecular gas tracers (Figure |9|. 

The combination of multiple clouds and nearly perpen- 
dicular orbital motion explains why the methods pre- 
sented in this paper do not match the maser parallax 
distances for these two clouds close to the Galactic cen- 
ter. We shade these sources in red in Figure [7] and Fig- 
ure [S] and analyze them separately in Table [4j We caution 
that extinction distances should always be combined with 
spectra of a dense gas tracer to detect multiple clouds 
along the line of sight and thereby identify distance esti- 
mates which are potentially unreliable. 

6.1.2. C23.66-0.13 

This source shows a significant discrepancy between 
the kinematic distances and the maser parallax distance. 
Although it is less well characterized as a region of high- 
mass star formation, there are several lines of evidence 
in support of this identification. The source was targeted 
for a methanol maser search based on its Infrared Astro- 
nomical Satellite (IRAS) colors; these observat ions re- 
vealed a rin g of 6.7 GHz methanol masers (Bartkiewicz 



et al.||2005| ). The discovery of additio nal objects with 
similar maser morphologies led Bartkiewicz et al. (2009) 
to propose a new class of "ring-like" methanol masers. 
The so urce is included in the Red MSX Source (RMS) 
Survey Urquhart et al. ( 2011, ) where it shows an spectral 
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Fig. 8. — Comparison of different distance estimates to the maser clouds. This is the same as Figure [T] except sorted by increasing maser 
parahax distance. Red shaded clouds show the two sources near the Galactic center where kinematic distances are unreliable and there is 
evidence for multiple clouds along the line of sight. The non-parallax methods (both kinematic and extinction) appear to be systematically 
biased toward deriving larger distances for near clouds {D < 2.5 kpc). 
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Fig. 9. — Multiple velocity components in dense gas tracers toward the two clouds closest to the Galactic center (right: G05. 88— 00.39, 
left: G09. 62+00. 20) indicat e the presence of multiple clouds along the line of sight, confounding extinction distance determinations here. 
CO (1-0) spectra are from jDame et al."2001) and show an abundance of velocity features toward these positions. The HCO+ (1-0) and 
H^^CO+ (1-0) spectra come from the MALT90 Survey (Foster et al. 2011). The high critical density (>10^ cm~^) of these molecules mean 
they are picking out just dense clouds. Vertical line shows the velocity used for kinematic distances. Both spectra show multiple velocity 
components in HCO+ (1-0). The optically thin H-'^'^CO+ (1-0) spectrum shows that the multiple velocity features seen in the HCO+ (1-0) 
spectra are not due to self-absorption. 
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energy distribution consistent with an HIT region and 
has detections of NH3 (1,1), (2,2) and (3,3) with Vlsr 
= 80.45 km s~^ toward the position. 

We detect a smah dark cloud associated with this 
maser position. Figure 10 shows the 0.12 mJy beam~^ 
contours of the BGPS data at this cloud overlaid on 
the 2MASS image. The ^^CO (1-0) spectrum from the 
GRS of this position shows a single strong line at 80 
km s~^ consistent with both the NH3 velocity and the 
central maser velocity and spatially consistent with the 
dark cloud. Because the dark cloud toward this source 
is so small, extinction distances are difficult. Only the 
UKIDSS data contains enough blue stars for a (highly 
uncertain) estimate of the distance from the Blue Num- 
ber Count method. The Red Giant Method does identify 
a cloud at 3.2^0.4 kpc by choosing a very small extrac- 
tion window around the maser position. However, this 
extinction feature is only 1 magnitude in Ay, far below 
our normal column density threshold for identifying a 
cloud; for this reason we do not consider this distance a 
reliable estimate. 

In conclusion, although G23.66— 0.13 is not as well- 
studied a region of high-mass star formation as the other 
objects considered in this study, it does appear to be as- 
sociated with a small dense cloud with a velocity consis- 
tent with its 6.7 GHz methanol masers, lending credence 
to the identification of this object as a young high-mass 
star. The other possible explanation is that the maser 
parallax distance is measuring the distance to an asymp- 
totic giant branch (AGB) star which is spatially coinci- 
dent with the dark cloud. AGB stars can power masers, 
though there is no known case of an AGB star power- 
ing a 6.7 GHz methanol maser. We conclude that the 
AGB explanation is unlikely and the significant differ- 
ence between the kinematic distance and maser parallax 
distance remains unexplained. 

6.1.3. G34-39+0.22 

For the source G34.39+0.22, the extinction distances 
presented herein have small uncertainties and are con- 
sistent with the kinematic distances for this cloud, but 
are inconsistent with the maser parallax distance (the 
2MASS-based Blue Number Count method is consistent 
with the maser parallax distance due solely to its large 
uncertainty) . The Red Giant extinction distance method 
identifies two clouds towards this line of sight, but the 
first cloud is a rather tenuous identification and spec- 
tra of dense gas tracers toward thi s source (e.g CS (2-1), 
HCO+ (1-0); Sanhueza et al. 12010 ) do not show a second 
velocity component. Using either of the rotation curves 
in this work, the cloud would need to have Ylsr "^20 km 
s~^, a discrepancy of 37 km s~^, in order for the kine- 
matic distance to match the maser parallax distance. 

There is weak CO emission around 27 km s~^ seen 
in both the GRS ^^CO (1-0) spectra and Atacama 
Pathfinder Experiment (APEX) ^^CO (3-2) spectra 
(Sanhueza, private communication), but it seems highly 
unlikely that t h e mul tiple maser spots measured in Ku- 
rayama et al. ( 2011| ) arise from this weak foreground 
cloud. 'I'he main cloud is quite large and well studied; 
the masers lie in known regions of high-mass star forma- 
tion which lie precisely along the dark filament we iden- 
tify for use in the extinction distance. The near-infrared 



data does not appear abnormal in any way. 

This source remains a puzzle. It is possible that the 
maser parallax distance is incorrect, or at least that 
the errors are significantly underestimated. Because this 
source is very far south for the VERA array (+1 ° of dec- 
lination), the declination does not constrain t he parallax, 
which come s only from fitting right ascension ( Kurayaina] 



et al. 



2011). In addition, the model for the right ascen- 



sion data (annual parallax plus a linear proper motion) 
is a poor fit to the observations, with 8 out of 21 observa- 
tions excluded from the fit for the three maser features 
used. Additional maser parallax measurements of this 
source would be highly desirable. 

7. SUMMARY 

New extinction distance measurements have been 
made toward 11 dark molecular clouds for which di- 
rect distances are available from maser parallax measure- 
ments. The two extinction methods are the Blue Number 
Count method and the Red Giant method. In addition, 
we have determined kinematic distances to these clouds 
using dense gas tracers and two different rotation curves 
(Clemens 1985; Reid et al. 2009). 

'I'he Blue Number Count method presented here uses 
auxiliary information to identify the boundaries of a 
cloud within which there is enough column density to 
cleanly separate all foreground stars from background 
stars. The number of blue foreground stars within these 
boundaries is then compared to a Galactic model of the 
stellar distribution to estimate the distance to the cloud. 
Because the surface density of blue stars must be esti- 
mated accurately we have used injection and recovery of 
synthetic stars to estimate the completeness as a func- 
tion of magnitude for each region (for both 2MASS and 
UKIDSS). 

The Blue Number Count method requires some addi- 
tional information in order to determine if there are mul- 
tiple clouds along the line of sight; if there are multiple 
clouds along the line of sight the method can only esti- 
mate a distance to the nearest cloud. The Blue Number 
Count method also only works for clouds with sufficiently 
high column density over a large enough region to have a 
significant number of blue foreground stars in whichever 
near-infrared catalog is being used. Possible sources of 
systematic error include multiple clouds along the line of 
sight, inhomogeneities in the stellar distribution, errors 
introduced by the completeness estimation and photome- 
try routines, and deficiencies in the Galactic model used. 

Our reliance on stellar number density means that the 
Blue Number Count met hod te sts a different aspect of 
the Besancon model than Marshall et al.| ( |2009^ ), which 
uses only color information and discards number density 
information by normalizing within each color bin. Mar-l 



shall et al. ( 2009J found a significant (1.5 kpc) systematic 
offset between extinction distances and kinematic dis- 
tances in the fourth quadrant. Although this work does 
not directly address this offset, since none of the clouds 
with maser parallax distance determinations are in the 
fourth quadrant, the use of complimentary information 
from the Besancon m odel will allow us to examine the 
Marshall et al. (2009) result in the fourth quadrant in 
future work. 

The Red Giant method used here is similar to that 



presented in Stead fc Hoare^ (2010} , whereby individual 
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Fig. 10. [Left] 2MASS three-color image (J = blue, H = green, K = red) of G23. 66-0.13. BGPS contours at 0.12 mjy beam-^ are 
in yellow, and the position of the maser source is shown by the cyan cross. [Right] GRS -"^"^CO (1-0) spectrum at the position of the maser 
source G23.66-0.13. The strongest emission is at +80 km s"-*^, consistent with the RMS NH3 velocity for this source ( Urquhart et al.|201l| 
and the central velocity of the methanol masers in this source (+82 km s"-*^), confirming the association of the maser source with the cloud 
identified in BGPS and GRS. 



extinction measurements are made to color-selected gi- 
ant stars along the same line of sight as molecular clouds. 
Following this, extinction histograms are constructed and 
molecular clouds can be identified as large gaps in what 
would otherwise be a relatively smooth distribution. Us- 
ing an existing extinction-distance relation (EDR) the 
extinction measurements can be converted to reliable dis- 
tances. However this method suffers at near distances 
where an inadequate number of giant stars in front of 
the cloud are available. The EDRs used here are those of 



IMarshall et al. (2006). As they gridded the Galaxy into 
15'xl5' tiles to derive their EDRs, the spatially closest 
tile to each cloud has been used. The reliability of this 
method could therefore be improved if line of sight spe- 
cific EDRs were derived using the very same stars used 
to build the extinction histogram. 

In some cases multiple molecular clouds have been 
identified along the same line of sight; such cases pose a 
challenge for extinction distance determinations. In this 
sample, multiple molecular clouds are typically observed 
towards the Galactic center where the surface density of 
clouds is high. The Blue Number Count method can 
only determine the distance to the first cloud along a 
line of sight, and may fail even for the foreground cloud 
if the multiple clouds along the line of sight cause us to 
misidentify the region which is dense enough to produce 
a clear separation between background and foreground 
stars. The Red Giant method can, with sufficient num- 
ber of stars, identify multiple clouds along the line of 
sight, a clear advantage for this method in crowded re- 
gions. 

Our individual distance determinations are presented 
in Table [3j The extinction distance methods agree to 
within 2(7 with the maser parallax distances between 66% 



and 100% of the time (depending on method and input 
survey), and between 85% and 100% of the time outside 
of the crowded Galactic center (Table |4|. Extinction 
distance estimates provide better agreement with maser 
parallax distances than the kinematic distances obtained 
using either of two different rotation curves. Extinc- 
tion distance methods provide an important, indepen- 
dent, estimate of the distance to dark molecular clouds. 
These methods do fail occasionally, but normally there is 
a good explanation such as multiple clouds along the line 
of sight, which can be established through auxiliary in- 
formation. In a small number of cases it is possible that 
that maser parallax distance is incorrect, or is reported 
with error bars which are too small. In the absence of 
a maser parallax distance, the combination of both ex- 
tinction and kinematic distances will provide the most 
secure measurement of the distance to a cloud. 
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Fig. A1. — Similar to Figure[T][2] andls] For each source, we show (left) the 2MASS three-color image (J=blue, if=green, K=red) with 
the cloud contour and the position of tne maser overlain. In the center we show the histogram of UKIDSS (where available, otherwise 
2MASS) J-K colors inside the cloud boundaries, and show the cutout at J-K = 1.5. On the right we show the UKIDSS (where available, 
otherwise 2MASS) completeness-corrected blue stellar density and how this compares to the Galactic model. 



Maser Extinction Distance Comparison 



17 




18 



Foster, Stead, Benjamin, Hoare & Jackson 




Fig. A3. — Same as Figure [AT] 
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Fig. B1. — Similar to Figurel6](a) and (b), each A^^ histogram has been labeled with the name of the target cloud - the 2MASS 
histograms have been denotedwith a '2M'. The near and far sides of each molecular cloud are marked with dashed lines, derived 
from the la deviations of the skewed Gaussians fitted to each histogram (solid curve). Distance- Ay plots have also been created, 
similar to Fig. [g] (c), showing the spatially closest EDR to each cloud (solid thin line). The two sets of dashed lines illustrate 
the near and far sides of a molecular cloud. 
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Fig. B3. Same as Fig. IbT 



